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We study the dynamical properties of the Mermin model, 
a simple quantum dissipative model with a monochromatic 
environment, using analytical and numerical methods. Our 
numerical results show that the model exhibits a second order 
phase transition to a localized state before which the system 
is effectively decoupled from the environment. In contrast 
to the spin-boson model, the Mermin model exhibits an "or- 
thogonality catastrophe," denning the critical point, before 
dissipation has destroyed all coherent behavior. An analytic 
approach based on the Liouvillian technique, though success- 
ful in describing the phase diagram of spin-boson and related 
models, fails to capture this essential feature of the Mermin 
model. 

PACS numbers: 03.65. Bz, 42.50.Lc 



I. INTRODUCTION 

The formal analogies between phase transitions in 
one-dimensional systems, decoherence in the spin-boson 
Hamiltonian (SBH), and screening in Kondo systems 
have been important in sorting out the subtle behavior of 
dissipative quantum systems. In this paper we undertake 
an analytic and numerical study of the Mermin model, a 
model closely related to the spin-boson model, but pos- 
sessing a monochromatic environmental spin bath. 

The spin-boson model is a universal model for a two 
level system interacting weakly with an environmen- 
tal bath, which has a continuous spectrum Al- 
though both monochromatic and polychromatic baths 
bring about a localization transition in the two level sys- 
tem, we point out, in this paper, an important difference 
between the dynamics of the two level system interacting 
weakly with these two kinds of baths. If the dynamics 
of the two level system is characterized by a (pseudo) 
spin susceptibility, increasing the coupling to the environ- 
ment broadens the resonance and shifts the resonant fre- 
quency somewhat. When the environment has a continu- 
ous and ohmic spectrum (the SBH) , the system spin hrst 
loses all coherent behavior at some intermediate coupling 
strength. The dynamics of the two level system becomes 
overdamped at this point, but the localization transition 
does not occur until a stronger coupling is reached. The 



critical coupling at which the two level system is ulti- 
mately localized may be identified with an "orthogonal- 
ity catastrophe." Thus, before the environmental overlap 
effects become catastrophic, they are finite and lead to a 
continuous softening of the system energy scale, as well 
as a complete inelastic broadening to an incoherent state 
(the so called "Toulouse point"). 

On the other hand, when the spectrum of the envi- 
ronment is monochromatic (the Mermin model), the re- 
sulting dynamics are different. As the coupling to the 
environment increases, the only effect is to broaden the 
coherent resonance. At a critical coupling, a redistribu- 
tion of spectral weight is initiated and spectral weight 
flows from the coherent resonance to a localized feature 
at zero frequency. For coupling strengths beyond this 
value, spectral weight is increasingly transferred from the 
resonant peak to the localized feature. Thus, unlike the 
SBH, the Mermin model does not show any overdamped 
feature before the localization transition sets in. When 
localization does set in (at the critical point), the coher- 
ent peak continues to exist, albeit with reduced spectral 
weight. This is our main result. 

The paper is organized as follows. In section ||, we 
introduce the Mermin model and study the dynamics of 
the model using an analytical technique. This technique 
(based on Liouvillian operator methods) has been applied 
successfully in the past, to a variety of dissipative mod- 
els, including the SBH. In section III, we study the model 
numerically, and show that the model exhibits a localiza- 
tion transition at a critical value of the coupling strength. 
The final section contains some concluding remarks, sum- 
marizing our results and suggesting future directions of 
study. 



II. LIOUVILLIAN DYNAMICS OF THE MERMIN 
MODEL 

In this section, we introduce the Mermin model and 
study the dynamics of the model analytically. The model 
introduced by Mermin is obtained from the well known 
spin-boson Hamiltonian, by retaining only the lowest two 
levels of each harmonic oscillator comprising the environ- 
mental bath 101. The Hamiltonian therefore describes a 
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"system" spin-1/2, represented by the Pauli matrices a, 
coupled to a set of N "environment" spin-1/2 degrees 
of freedom, {s.,}. The resulting Hamiltonian, called the 
Mermin model is given by 
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In this paper, we investigate a large N variant of the 
model (0). As first noted by Mermin, the model can be 
solved approximately in the N — > oo limit, to demon- 
strate the correspondence between the onset of localiza- 
tion of the system spin and a second order phase tran- 
sition [Q]. In this large N model, the environment spins 
are summed to one big O(N) spin, S, and contributions 
from the sectors of Hilbert space with total spin S < N/2 
are ignored. The Hamiltonian then becomes 
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As Mermin points out, the advantage of the Hamil- 
tonian (||) is that in the limit N — > oo, the environ- 
ment spins may be replaced in the Hamiltonian by the x 
and z components of a classical spin angular momentum: 
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sulting Hamiltonian may be diagonalized and its ground 
state eigenvalue, E (9), given by 
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minimized with respect to 9. 

The critical behavior of the model may now be seen 
by examining the ground state energy Eq, as a function 
of A. The ground state energy, Eq, bifurcates at a finite 
value of A going from a singlet, non-degenerate root for 
A > X 2 /2uj to a doubly degenerate set for A < A 2 /2w: 
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In the former case, the environment is decoupled from the 
system and always points along the z-axis to minimize its 
"Zeeman" energy corresponding to the last term of (^). 
In the latter case, the environment and the system spins 
are frozen in two distinct orientations with degenerate 
energies. The ground state wave function is still in a 
superposition a\+) + /3|— ), but is no longer an evenly 
weighted one (a = /3); rather, the two roots correspond 
to the system predominantly in |+) or |— ). 

Despite producing a localization transition, this mean 
field solution is necessarily incomplete, for, it ignores 
quantum fluctuations of the bath; i.e., it provides no 
mechanism for dissipation. As a first step in this direc- 
tion, we study the dynamics of the bath using the Liou- 
villian operator method. This method has been applied 



successfully to study a wide variety of dissipative sys- 
tems. Using the Liouvillian approach, Shao and Hanggi 
H studied a spin-spin bath system similar to the Hamil- 
tonian (0) above, but with an ohmic spectrum and found 
a phase diagram virtually identical to that of the SBH. 
Dattagupta et al Q, studying the SBH, also employed 
a Liouvillian approach and found close agreement with 
other analytic approaches to the SBH jl]|5| . 

To try to gain some insight into the Mermin model be- 
yond the mean field level, we have adapted the resolvent 
operator calculation of Shao and Hanggi || on the spin- 
spin bath model. In what follows, we just outline the 
basic steps in the calculation and describe our results. A 
detailed review of the resolvent operator formalism can 
be found elsewhere pL 

The Hamiltonian (|2j) is rotated by R = exp — ia z 9^S y 
to diagonalize the environment spin for the two orienta- 
tions of the system spin. We then get 
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To study the dynamics of the system spin, a suitable 
quantity to calculate is (a z (t)} where (. . .) denotes quan- 
tum average at zero temperature. For a state prepared 
initially at t — 0, with the system and environment spin 
fully polarized, the time dependent average of the system 
spin is given by 

P (t) = {S T \R^ z (t)R\ T S) , 

where the first index of the ket is the system state and 
the second is the environment. The spin operator a z (t) is 
in the Heisenberg representation, and its time evolution 
is governed by the transformed Hamiltonian, H'. We 
now introduce the Liouvillian time evolution operator L, 
which is defined generically for a Hamiltonian H, Hilbert 
space {\j)} and operator O as follows: 

(i\O(t)\j)=J2(ij\e iU \i , f){i'\O(0)\f) ■ 

i',3' 

The time evolution operator Q(t) = exp (iLt) is a fourth 
rank operator whose state space vectors are denoted by 
a ket-parenthesis, All dynamical information is con- 
tained in Q and the Laplace transform of Q 

(jfikv\Q(z)\l[/mi/) = (j[ikv\ \lfi'mi/') , 

z — iL 

is the most convenient form to calculate. For the initial 
conditions specified above, it follows that the Laplace 
transform of pit) is given by: 

p(z) = ^2 aja*k{ k tit \Q{z)\mfimii)n (5) 

mfj>jk 
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where the a's are the rotation group matrix elements, 



(j T \R\ T S) 



Following Dattagupta et al. Q , we concentrate on the 
"environment-averaged" time evolution operator defined 
by performing the environment sums in ([5]) 

(fj,v\{Q(z))\(i! V) = ^a J a* k {k^jv\Q{z)\m^mv l ) 

kjm 

The Liouvillian is now split into perturbation (Li) and 
environment (L e ) parts, corresponding to the two terms 
in equation (^). Then the environment averaged time 
development operator may be computed perturbatively 
as 



(H(Q(*V) = GH 
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(7) 
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and G is the Liouvillian propagator for the environment 
degrees of freedom, G~ x — z — iL e . 

We can now calculate the above terms. The results 
through second order in A for the three self-energies in 
the order appearing in equation (Q) (notated Mi, M21, 
M22) are given below: 



Mi = iA(S\R 2 \S}(5 ® a x - a x ® 8) 



(8) 



Fig. ([j]). Since — e„ depends only upon k—n, n = k+r, 
may be written in terms of a relative index r and the k 
sum performed. For large S, the matrix elements, as 
functions of r, are peaked at the most probable values of 
r. Now, the two self energies may be written in the large 
S, continuum limit as follows: 
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FIG. 1. Fluctuation matrix elements depicted for processes 
contributing to B(z) and C(z). 

These self energies cannot be calculated in closed form; 
however, they can be evaluated in several limits that are 
physically significant. For C(z) we find, 



M 2 i = ~B{z)8 ®{8 + o z ) - -^B*{z){8 + a z )®8 
~C[z)6 ®{8- <j z ) - ^C*{z){8 ~a z )®8 



+ 1 -{B(z) + B*{z))S- 



+ -(C(z) + C*(z))S + ®S+ 



M 22 = l*-\(S\n! i \S)\ 2 (8®5-o*®o*) 



(9) 
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The outer product notation a <g> [3 corresponds to 
aififi' fiw 1 ■ The z dependences of M21 are given by 

kn 
kn 

where ekn = £fc — e« is the bare energy difference. The 
fluctuation matrix elements are depicted graphically in 



ReC(z) 
ReC(z) 
and for B(z), 

ReB(z) 



A 2 

— z << A 

A 

A 2 

z 



>> A or A << 



A 



ReB(z) ~ A 2 
A 2 

RcB(z) 

z 



z 2 + e 2 



z « A 
A << oj 

z » A 



where e = uj (cos 8 — cos3#o)- 

To compute the Laplace transform of p(t) the self en- 
ergies (||), (|^), and ( ^Tj[ ) must be inserted in eqn. (^) 
and the fourth rank operator inverted. The general form 
is complicated, but for small z, the analytic structure 
reduces to 



Viz) 



ReC(z) + z 



z(z + ReC(z) +ReB(z)) 
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In the limit of large environmental coupling, the real pole 
dominates and the behavior for p(t) is exponential decay; 
at zero temperature, the system prepared in one state will 
leak into the other state. Thus, we get, 



Pit) 
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For small environmental coupling, the pole in p(z) is 
imaginary and p{t) oscillates, given by 



p(t)= cos At A << A 



(12) 



The results Eq. |fl]) and (|T|) summarize the two ex- 
tremal behaviors of p(t), but give no information about 
the critical point between. 



III. NUMERICAL RESULTS FOR THE MERMIN 
MODEL 

There are two main differences between the predictions 
of the Liouvillian perturbative scheme and the exact nu- 
merical behavior. First, the Liouvillian formalism does 
not capture the phase transition observed in finite size 
scaling of the numerical results (as well as in the mean 
field solution) described below. This result is interesting, 
given that the Liouvillian formalism does describe the 
phase transition at a = 1 of the spin-boson [Q and spin- 
spin bath H models. Consequently, one would expect it 
to work for the simpler monochromatic spectrum of the 
Mermin model. Secondly, the Liouvillian formalism pre- 
dicts exponential decay in p(t) ~ e~~ tA / A for large envi- 
ronmental coupling A. The corresponding spectral func- 
tion is Lorentzian centered at v = 0, similar to that of the 
spin-boson model at the Toulouse point a — 1/2. How- 
ever, in our numerical results we will show that the spec- 
tral function has only a finite energy (~ 0(A)) feature 
and an exponentially small frequency (~ 0{Ae~ N / 2 )) 
delta function feature (for environmental coupling larger 
than the critical value.) 

Let us now summarize briefly, our numerical findings. 
The transition to decoherence is induced by strong cou- 
pling to environmental degrees of freedom (large A), or, 
equivalently, a small system energy scale (small A). Be- 
low a critical coupling, A c , all spectral weight of the dy- 
namical susceptibility 



X » = Im 



dt(0\[a(t),a(0)}\0}e(t)e l 



(13) 



resides in the principal resonance of the two level sys- 
tem at an energy A. When the coupling is larger than 
A c , a new exponentially small energy scale 0(Ae~ N / 2 ) 
emerges, associated with a broken symmetry, (a) =/= 0, 
in the thermodynamic limit N — ► oo. This feature 
corresponds to tunneling modified by a Franck-Condon 



type overlap factor. As the coupling is increased, spec- 
tral weight is shifted continuously to the "near-zero" fre- 
quency channel. The weight of the delta- function, 5(0 + ), 
is simply the order parameter, |(0|er|0)| 2 (or, more ex- 
actly, |(l|cr|0)| 2 , for large but finite N). However, the dy- 
namical susceptibility obeys a sum rule, implying that an 
incompletely formed broken symmetry state leaves some 
spectral weight at the position of the principal resonance, 
set by A. 

To calculate the dynamical susceptibility we consider 
the response of the system (0) to an external field h(t) = 
hcos vt: 



H, 



ext 



-/i— cos VI 



(14) 



The Hamiltonian (g) was diagonalized and the com- 
plete order parameter matrix, cry = (i|ct|j), was com- 
puted. It is natural to work with the dimensionless cou- 
pling constant, k = 2luA/X 2 ; The critical coupling im- 
plied in the N — > oo calculation is then k = k c = 1. 




FIG. 2. Evolution of spectral weight in the Mermin Model. 
x"(^) is shown for four values of k decreasing below the crit- 
ical value, k c . Despite the change in the position of the prin- 
cipal resonance (due to changing A), spectral weight emerges 
at v ~ abruptly at n c and the two features remain distinct. 

Fig. (|J) shows the behavior of x"(f) as the system 
energy scale A is reduced, corresponding to a range 
k = 0.6« c — l.2n c . These computations were performed 
for N — 80. Computations on larger spin environments 
(up to N = 300) suggest that there is little change be- 
yond N = 80. An artificial broadening (5 = 0.01) was 
added to make the features more visible. As seen in this 
sequence, n = n c is marked by the appearance of the 
Franck-Condon resonance at an exponentially small en- 
ergy scale (v\q ~ 10~ 3 ). The spectral weight remaining 
at 0(t) when k > n c clearly exhibits inelastic broaden- 
ing, although the resonance essentially disappears before 
becoming critically damped. 
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The counterpart, at finite N, to (a) ^ at infinite N 
is the spectral weight at the Franck-Condon resonance, 
oxo 7^ 0. Fig. 3 demonstrates that uio behaves as ex- 
pected in a second order phase transition; the slope at 
n = grows with increasing N. The matrix, try, and 
consequently, x"(f), obey a sum rule: 



(15) 



An incompletely developed broken symmetry ( | cr 10 1 2 < 
1) means that spectral weight remains at the principal 
resonance at 0(A). 
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FIG. 3. The order parameter, do, is plotted as a func- 
tion of re and exhibits second order phase transition behavior, 
re was varied by changing A. Comparison of N = 80 and 
N — 300 data show the slope steepening at the critical point. 
When re > re c , there is no weight in the Franck-Condon reso- 
nance other than that attributed to the rounding of the phase 
transition at finite N. 

The quantum fluctuations missed by the mean field 
analysis but captured in the numerics may be qualita- 
tively divided into two types: overlap effects and inelas- 
tic effects. In the former type, the system in the |+) 
state leaves an "imprint" upon the environment which, 
being nearly orthogonal to the imprint left by system 
state |— ), reduces the tunneling amplitude by an expo- 
nentially small factor. The new tunnel splitting is thus 
reduced by a Franck-Condon type overlap factor: 



A* = A(0oK(20o)|0o) 



A cos^ 0o 



1 



'(l + A 2 /w 2 ) w / 2 
(16) 

where R y (20o) is the rotation operator. The reso- 
nance then shifts from O(A) to a smaller energy scale 
0(Ae _JV / 2 ). In the inelastic type fluctuation, the sys- 
tem forces the environment to make transitions to an 
excited state dissipating an amount energy A 2 /u> per pe- 
riod. The resonance is then broadened by 0(A 2 /w) but 
remains nominally at O(A). 



These results show that the spin dynamics of the Mer- 
min model is significantly different from that of the SBH 
with an ohmic bath. The former, as our results show, 
only exhibits overlap effects at the critical coupling, but 
not before. In the SBH, the coupling to the environ- 
ment produces finite overlap effects for all couplings 
< a < 1/2, and a Franck-Condon type reduction in 
energy scale, given by: 



A cff = A(— )i 



(17) 



At the same time, inelastic effects decohere the resonance 
leading to the Toulouse limit (at a = 1/2), prior to lo- 
calization (which occurs at a = 1) , and corresponding to 
a complete inelastic broadening of the resonance. 
In the Mermin model, the quantum resonance of the sys- 
tem spin, although damped, remains intact through the 
localization transition. 



IV. CONCLUSION 

In this paper we have compared the dynamical be- 
havior of a quantum dissipative system possessing a 
monochromatic spectrum (the Mermin model) with that 
of the spin-boson model which possesses a continuous 
bath spectrum. We find an important difference be- 
tween the dynamical behavior of the two models, with 
the Mermin model reaching a localized state, as environ- 
mental coupling is increased, without the primary res- 
onance becoming overdamped. When the critical cou- 
pling is reached, spectral weight is tranferred across a 
gap of O(A) to a zero frequency delta function, indicat- 
ing the presence of a localized state. In effect, the Mermin 
model exhibits no finite environmental overlap effects 
until an infinite one — the orthogonality catastrophe — 
localizes the system. In contrast, the primary resonant 
feature in spin-boson model must soften to zero frequency 
and become completely overdamped (the Toulouse point 
a = 1/2) before a localized spectral feature appears at 
a = 1. Furthermore, the Liouvillian technique, which 
has been successful in capturing the localization phase 
transition in the SBH and related models fails for the 
Mermin model. 

The Mermin model is a prototype system for a weakly 
coupled 0(1/ N) monochromatic spin bath, and judging 
from the results of Hanggi, a weakly coupled monchro- 
matic oscillator bath is likely to exhibit the same behav- 
ior. In contrast, the spin-boson model employs a feature- 
less, power-law spectrum of envrionmcntal oscillators and 
exhibits universal behavior that depends only upon the 
exponent of the spectrum. This behavior, in the ohmic 
case, includes a localization transition, but one with a 
qualitatively different character than that of the Mermin 
model. It is natural, then, to ask under what condi- 
tions does a quantum environment become sufficiently 
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"featureless" to fit into the universal framework of the 
spin-boson model? 
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